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Abstract 

We investigate the global texture model of structure forma- 
tion in cosmogonies with non-zero cosmological constant for 
different values of the Ifubble parameter. We find that the 
absence of significant acoustic peaks and little power on large 
scales are robust predictions of these models. However, from a 
careful comparison with data we conclude that at present we 
cannot safely reject the model on the grounds of present CMB 
data. Exclusion by means of galaxy correlation data, requires 
assumptions on biasing and statistics. New, very stringent 
constraints come from peculiar velocities. 

Investigating the large- TV limit, we argue that our main 
conclusions apply to all global 0{N) models of structure for- 
mation. 



I. INTRODUCTION 

Recently, a lot of effort has gone into the determina- 
tion of cosmological parameters from measurements of 
cosmic microwave background (CMB) anisotropics, es- 
pecially in view of the two planned satellite experiments 
MAP and PLANCK Q. However, we believe it is im- 
portant to be aware of the heavy modeling which enters 
these results. In general, simple power law initial spectra 
for scalar and tensor perturbations and vanishing vector 
perturbations are assumed, as predicted from inflation. 
To reproduce observational data, the composition of the 
dark matter and the cosmological parameters as well as 
the input spectrum and the scalar to tensor ratio are 
varied (|]. 

We want to take a different approach: We modify the 
model for structure formation. We assume that cosmic 
structure was induced by scaling seeds. Using a simplified 
(and not very accurate) treatment for the photon prop- 
agation, we have already shown that some key observa- 
tions can be reproduced within a very restricted family of 
scaling seed models ^ . Here we want to outline in detail 
a more accurate computation with a fully gauge- invariant 
Boltzmann code especially adapted to treat models with 
sources. In this paper we follow the philosophy of a gen- 
eral analysis of scaling seed models motivated in Ref. . 

Seeds are an inhomogeneously distributed form of mat- 
ter (like e.g. topological defects) which interacts with the 



cosmic fluid only gravitationally and which represents al- 
ways a small fraction of the total energy of the universe. 
They induce geometrical perturbations, but their influ- 
ence on the evolution of the background universe can be 
neglected. Furthermore, in first order perturbation the- 
ory, seeds evolve according to the unperturbed spacetime 
geometry. 

Here, we mainly investigate the models of structure for- 
mation with global texture. This models (for ^matter = 
1) show discrepancies with the observed intermediate 
scale CMB anisotropics and with the galaxy power spec- 
trum on large scales Recently it has been argued 
that the addition of a cosmological constant leads to bet- 
ter agreement with data for the cosmic string model of 
structure formation . We analyze this question for the 
texture model, by using ab initio simulation of cosmic 
texture as described in Ref. We determine the CMB 
anisotropics, the dark matter power spectrum and the 
bulk velocities for these models. We also compare our 
results with the large-A^ limit of global 0{N) models, 
and we discuss briefly which type of parameter changes 
in the 2-point functions of the seeds may lead to better 
agreement with data. 

We flnd that the absence of significant acoustic peaks 
in the CMB anisotropy spectrum is a robust result for 
global texture as well as for the large-iV limit for all 
choices of cosmological parameters investigated. Further- 
more, the dark matter power spectrum on large scales, 
A ^ 20/i~^Mpc, is substantially lower than the measured 
galaxy power spectrum. 

However, comparing our CMB anisotropy spectra with 
present data, we cannot safely reject the model. On large 
angular scales, the CMB spectrum is in quite good agree- 
ment with the COBE data set, while on smaller scales we 
find a significant disagreement only with the Saskatoon 
experiment. Furthermore, for non-satellite experiments 
foreground contamination remains a serious problem due 
to the limited sky and frequency coverage. 

The dark matter power spectra are clearly too low on 
large scales, but in view of the unresolved biasing prob- 
lem, we feel reluctant to rule out the models on these 
grounds. A much clearer rejection may come from the 
bulk velocity on large scales. Our prediction is by a fac- 
tor 3 to 5 lower than the POTENT result on large scales. 

Since global texture and the large- TV limit lead to very 
similar results, we conclude that all global 0{N) models 
of structure formation for the cosmogonies investigated 
in this work are ruled out if the bulk velocity on scales of 
50/i^^Mpc is around 300km/s or if the CMB primordial 
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anisotropics power spectrum really shows a structure of 
peaks on sub-degree angular scales. 

This paper is the first of a series of analyses of mod- 
els with scaling seeds. We therefore fully present the 
formalism used for our calculations in the next section. 
There, we also explain in detail the eigenvector expan- 
sion which allows to calculate the CMB anisotropics and 
matter power spectra in models with seeds from the two 
point functions of the seeds alone. This section can be 
skipped if the reader is mainly interested in the results. 
Section 3 is devoted to a brief description of the numeri- 
cal simulations. In Section 4 we analyze our results and 
in Section 5 we draw some conclusions. Two appendices 
are devoted to detailed definitions of the perturbation 
variables and to some technical derivations. 

Notation: We always work in a spatially flat Fried- 
mann universe. The metric is given by 

ds^ = a{tf{dt^ - Sijdx'dx^) , 

where t denotes conformal time. 

Greek indices denote spacetime coordinates (0 to 3) 
whereas Latin ones run from 1 to 3. Three dimensional 
vectors are denoted by bold face characters. 



II. THE FORMALISM 

Anisotropics in the CMB are small and can thus be 
described by first order cosmological perturbation theory 
which we apply throughout. We neglect the non- linear 
evolution of density fluctuations on smaller scales. Since 
models with seeds are genuinely non-Gaussian, the usual 
numerical N-body simulations which start from Gaussian 
initial conditions cannot be used to describe the evolution 
on smaller scales. 

Gauge-invariant perturbation equations for cosmolog- 
ical models with seeds have been derived in Refs. [^,^. 
Here we follow the notation and use the results presented 
in Ref. Definitions of all the gauge-invariant pertur- 
bation variables used here in terms of perturbations of 
the metric, the energy momentum tensor and the bright- 
ness are given in Appendix A for completeness. 

We consider a background universe with density pa- 
rameter fio = + = 1, consisting of photons, cold 
dark matter (CDM), baryons and neutrinos. At very 
early times z 3> Zrfec ^ 1100, photons and baryons form 
a perfectly coupled ideal fluid. As time evolves, and as 
the electron density drops due to recombination of pri- 
mordial helium and hydrogen, Compton scattering be- 
comes less frequent and higher moments in the photon 
distribution develop. This epoch has to be described by 
a Boltzmann equation. Long after recombination, free 
electrons are so sparse that the collision term can be ne- 
glected, and photons evolve according to the coUisionless 
Boltzmann or Liouville equation. During the epoch of 
interest here, neutrinos are always coUisionless and thus 
obey the Liouville equation. 



In the next subsection, we parameterize in a completely 
general way the degrees of freedom of the seed energy mo- 
mentum tensor. Subsection B is devoted to the pertur- 
bation of Einstein's equations and the fluid equations of 
motion. Next we treat the Boltzmann perturbation equa- 
tion. In Subsection D we explain how we determine the 
power spectra of CMB anisotropics, density fluctuation 
and peculiar velocities by means of the derived pertur- 
bation equations and the unequal time correlators of the 
seed energy momentum tensor, which are obtained by 
numerical simulations. In Subsection E we give the ini- 
tial conditions and a brief description of our Boltzmann 
code. 



A. The seed energy momentum tensor 

Since the energy momentum tensor of the seeds, 
Qfxv, has no homogeneous background contribution, it is 
gauge invariant by itself according to the Stewart- Walker 
Lemma ]To[ |. 

Ofj_i, can be calculated by solving the matter equa- 
tions for the seeds in the Friedmann background geometry 
(Since O^j, has no background component it satisfies the 
unperturbed "conservation" equations.). We decompose 
0^1^ into scalar, vector and tensor contributions. They 
decouple within linear perturbation theory and it is thus 
possible to write the equations for each of these contribu- 
tions separately. As always (unless noted otherwise), we 
work in Fourier space. We parameterize the scalar (S) 
vector (V) and tensor (T) contributions to 9^i/ in the 
form 



(1) 
(2) 

(3) 

(4) 
(5) 
(6) 



Here M denotes a typical mass scale of the seeds. In the 
case of topological defects we set M = rj, where rj is the 
symmetry breaking scale The vectors u;^"^ and w^^^ 

are transverse and t^J'^ is a transverse traceless tensor, 





- M^fp 










ifp + ^k'^f7T)Sji - kjkif^ 




Jl f2 (f) 

= M w )j 






KkM'^'+hw^^') 









fcV. 







From the full energy momentum tensor O^^ which may 
contain scalar, vector and tensor contributions, the scalar 
parts /„ and /^r of a given Fourier mode are determined 
by 



ik^Q. 
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On the other hand fy and /^r are also determined in terms 
of fp and fp by energy and momentum conservation, 



fp + k\h. + -{fp + 3fp) = , 



a 3 



Once is known it is easy to extract M^wj"'' = 8 
ikjM'^fy. For we use 



(7) 
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B. Einstein's equations and the fluid equations 

1. scalar perturbations 

Scalar perturbations of the geometry have two degrees 
of freedom which can be cast in terms of the gauge- 
invariant Bardeen potentials, \1/ and <I> For New- 
tonian forms of matter = — $ is nothing else than the 
Newtonian gravitational potential. For matter with sig- 
nificant anisotropic stresses, \1/ and — $ differ. In geo- 
metrical terms, the former represents the lapse function 
of the zero-shear hyper-surfaces while the latter is a mea- 
sure of their 3-curvature In the presence of seeds, the 
Bardeen potentials are given by 



2/l#2,„(^) 



Again, can also be obtained in terms of 

means of momentum conservation, 



„(") 



by 
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ifc^wl") = 
2 ' 



(9) 



The geometry perturbations induced by the seeds are 
characterized by the Bardeen potentials, and '^s, for 
scalar perturbations, the potential for the shear of the ex- 
trinsic curvature, for vector perturbations and the 

(s) 

gravitational wave amplitude, H^j , for tensor perturba- 
tions. Detailed definitions of these variables and their 
geometrical interpretation are given in Ref. [D (see also 
Appendix A). Einstein's equations link the seed pertur- 
bations of the geometry to the energy momentum tensor 
of the seeds. Defining the dimensionless small parameter 



we obtain 
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Eqs. ( [TTI ) to (|lj) would determine the geometric pertur- 
bations if the cosmic fluid were perfectly unperturbed. 
In a realistic situation, however, we have to add the fluid 
perturbations in the geometry which are defined in the 
next subsection. Only the total geometrical perturba- 
tions are determined via Einstein's equations. In this 
sense, Eqs. dnl) to (|l|) should be regarded as defini- 
tions for $s , *s , S^'-) and H^^f . 

A description of the numerical calculation of the energy 
momentum tensor of the seeds for global texture is given 
in Section III. 



(15) 
(16) 



where the indices s,m refer to contributions from a source 
(the seed) and the cosmic fiuid respectively. The seed 
Bardeen potentials are given in Eqs. ([l^) and (|l2|). 

To describe the scalar perturbations of the energy mo- 
mentum tensor of a given matter component, we use the 
variables Dg, a gauge-invariant variable for density fluc- 
tuations, V, the potential of peculiar velocity fluctua- 
tions, and n, a potential for anisotropic stresses (which 
vanishes for CDM and baryons). A definition of these 
variables in terms of the components of the energy mo- 
mentum tensor of the fluids and the metric perturbations 
can be found in Refs. ||l^ or ^ and in Appendix A. 

Subscripts and superscripts ^, c, 6 or denote the ra- 
diation, CDM, baryon or neutrino fluids respectively. 

Einstein's equations yield the following relation for the 
matter part of the Bardeen potentials p3[ 
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(17) 
(18) 



Note the appearance of <I> = "I>s + on the r.h.s. 
of Eq. (|l7|). Using the decompositions ([l5| , |l6| ) we can 
solve for $ and ^ in terms of the fluid variables and the 
seeds. With the help of Friedmann's equation, Eqs. (pi 
and (|l^) can then be written in the form 
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(i) k'^+Ax^ + Zxc + 2,Xb + Ax^ 
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+ -k 
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+4X.K) + (- 
3 \a 



- ) k {xjU^ + x^Ui,) 



(19) 



(20) 



Here we have normalized the scale factor such that 
a — 1 today. The density parameters always repre- 
sent the values of the corresponding density parameter 
today (Here . stands for c , 7 , 6 or j,.). To avoid any 
confusion, we have introduced the variables x, for the 
time dependent density parameters, 



Qc.ba 



(21) 
(22) 



The fluid variables of photons and neutrinos are ob- 
tained by integrating the scalar brightness perturbations, 
which we denote by Ms{t, k, n) and Afs{t, k, n) respec- 
tively, over directions, n, 
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(k • n)Msdn 



(23) 
(24) 
(25) 
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■ Ji(k-nf -^k')Msdn (26) 

(27) 

^9'^ = ^ / -^sdn = 1^0 (28) 
3i 
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(30) 

(31) 
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A systematic definition of the modes crj and Vj is given 
in the next subsection. 

The equation of motion for CDM is given by energy 
and momentum conservation, 



D^^^ +kVc = , 



Vc 



-Wc^k^ . 

a, 



(33) 
(34) 



i)^ + -kv^ = , 

V, - fci^^) = ki^ - $) , 

n(b) = ^ (^) 
9 4^9 



(35) 

(36) 

(37) 
(38) 



The conservation equations for neutrinos are not very 
useful, since they involve anisotropic stresses and thus 
do not close. At the temperatures of interest to us, 
T ^ IMeV, neutrinos have to be evolved by means of 
the Liouville equation which we discuss in the next sub- 
section. 

Once the baryon contribution to the baryon-photon 
fluid becomes non-negligible, and the imperfect coupling 
of photons and baryons has to be taken into account (for 
a 1% accuracy of the results, the redshift corresponding 
to this epoch is around z ~ lO''), we evolve also the 
photons with a Boltzmann equation. The equation of 
motion for the baryons is then 



kVh 



Vb 



, 
k^ 



(39) 



SQb 



-"[V.-Vb]. (40) 



The last term in Eq. (40) represents the photon drag 



force induced by non-relativistic Compton scattering, 
is the Thomson cross section, and denotes the num- 
ber density of free electrons. At very early times, when 
axrie ^ 1/^, the 'Thomson drag' just forces Vb = V~ 



7' 



which together with Eqs. (35) and ( |39D implies (^ 

An interesting phenomenon often called 'compensa- 
tion' can be important on super horizon scales, fci <C 1. 
If we neglect anisotropic stresses of photons and neutri- 
nos and take into account that OjDg) — 0{ktV) and 
0{V) = 0{km) for kt < 1, Eqs. M) and {Mj lead to 



0($) = O ((fc<)2$, - 2eU) 



(41) 



Hence, if anisotropic stresses are relatively small, e/jr <C 
$8, the resulting gravitational potential on super horizon 
scales is much smaller than the one induced by the seeds 
alone. One must be very careful not to over interpret this 
'compensation' which is by no means related to causality, 
but is due to the initial condition Dg , V ^t-to 0. A 
thorough discussion of this issue is found in Refs. [p^[-p5[. 
As we shall see in the next section, for textures $s and 
e/^ are actually of the same order. Therefore Eq. ( |4l| ) 
does not lead to compensation, but it indicates that CMB 
anisotropics on very large scales (Sachs- Wolfe effect) are 
dominated by the amplitude of seed anisotropic stresses. 



During the very tight coupling regime, z ^ Zdec, we may 
neglect the baryon contribution in the energy momentum 
conservation of the baryon-photon plasma. We then have 



The quantities which we want to calculate and compare 
with observations arc the CDM density power spectrum 
and the peculiar velocity power spectrum today 
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P{k) = {\D(^Hk,to)\' 
and 

Pv{k) = {\Vc{k,to)f) 



(42) 
(43) 



Here (• • •) denotes an ensemble average over models. 
Note that even though Dg and V are gauge invariant 
quantities which do not agree with, e.g., the correspond- 
ing quantities in synchronous gauge, this difference is 
very small on subhorizon scales (of order 1/kt) and can 
thus be ignored. 

On subhorizon scales the seeds decay, and CDM per- 
turbations evolve freely. We then have, like in inflation- 
ary models, 



P4k) = H^nl;,^P{k)k- 



(44) 



(see next section) and transfer it onto the baryons via 
the photon drag force. 



a J 61 lb 



(48) 



The photon vorticity is given by an integral over the vec- 
tor type photon brightness perturbation, Mv, 



4:71 



(49) 



where the integral is over photon directions, n. In terms 
of the development presented in the next section for k 
pointing in z-direction, we obtain 



^(V) (V) (V) (V) 
"1,2 ' ■ "2,2 + "2,0 ' ^ 



(50) 



2. vector perturbations 

Vector perturbations of the geometry have two degrees 
of freedom which can be cast in a divergence free vector 
field. A gauge-invariant quantity describing vector per- 
turbations of the geometry is S, a vector potential for the 
shear tensor of the {t =const.} hyper-surfaces. Like for 
scalar perturbations, we split the contribution to S into a 
source term coming from the seeds given in the previous 
subsection, and a part due to the vector perturbations in 
the fluid, 



Equivalently, we have for neutrinos 
47r ,/ 



(51) 
(52) 



The vector equations of motion for photons and neutrinos 
are discussed in the next section. 



3. tensor perturbations 



The perturbation of Einstein's equation for is 

2 A 



k 



(45) Metric perturbations also have two tensorial degrees 

of freedom, gravity waves, which are represented by the 
9| two helicity states of a transverse traceless tensor (see 

Appendix A). As before, we split the geometry perturba- 
tion into a part induced by the seeds and a part due to 
the matter fluids. 



(46) 



Here a;, is the fluid vorticity which generates the vec- 
tor type shear of the equal time hyper-surfaces (see Ap- 
pendix A). By definition, vector perturbations are trans- 
verse, 



S • k = S„ • k = S,, • k = a;. • k = . 



(47) 



It is interesting to note that vector perturbations in 
the geometry do not induce any vector perturbations in 
the CDM (up to unphysical gauge modes), since no geo- 
metric terms enter the momentum conservation for CDM 
vorticity, 

a 

U3c + -uic = , 
a 

hence we may simply set u^c = 0. This is also the case 
for the tightly coupled baryon radiation plasma. But as 
soon as higher moments in the photon distribution build 
up, they feel the vector perturbations in the geometry 



H. 



(m) 



(53) 



The only matter perturbations which generate gravity 
waves are tensor type anisotropic stresses which are 
present in the photon and neutrino fluids. The pertur- 
bation of Einstein's equation yields 



Hlf+2(^\H^^;'^ + k^H^f = 



^) (x,h(j)h-x.h(;)). 



(54) 



The relation between the tensor brightness perturba- 



tions A4t, A/t and the tensor anisotropic stresses, H 



(7) 



and Hj-J-* is given by 



(55) 
(56) 
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In terms of the development presented in the next section 
for k pointing in z-direction, we have 
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(57) 
(58) 
(59) 
(60) 



We find that the effect of anisotropic stresses of photons 
and neutrinos is less than 1% in the final result, and 
hence we have neglected them. 



1. scalar perturbations 

We expand the brightness A^s(k, n, t) in the form: 
00 

A^5(n, k, t) = Y,{-iY{2i + l)af ) (i, k)P,(M), (64) 
£=0 

where Pi denotes the Legendre polynomial of order I 
and (Tg is the associated multipole moment. An analo- 
gous decomposition also applies to the amplitude of po- 
larization anisotropy, Ai^{n, k, t), and we denote the as- 

( S) 

sociated multipole moment by . 

The Boltzmann equation for scalar perturbations in 
the photon brightness and polarization is Pjl6| 



C. The Boltzmann equation 

When particle interactions are less frequent, the fluid 
approximation is not sufficient, and we have to describe 
the given particle species by a Boltzmann equation, in or- 
der to take into account phenomena like coUisional and 
directional dispersion. In the case of massless particles 
like massless neutrini or photons, the Boltzmann equa- 
tion can be integrated over energy, and we obtain an 
equation for the brightness perturbation which depends 
only on momentum directions ||^. As before, we split the 
brightness perturbation into a scalar, vector and tensor 
component, and we discuss the perturbation equation of 
each of them separately]^ 



M = Ms + Mv +M'j 
and 

TV^Ws +7VV +7Vt. 



(61) 
(62) 
(63) 



The functions M and J\f depend on the wave vector k, 
the photon (neutrino) direction n and conformal time t. 
Linear polarization of photons induced by Compton scat- 
tering is described by the variable A^^*^) (the Stokes pa- 
rameter Q) depending on the same variables. We choose 
for each k-mode a reference system with z-axis parallel 
to k. For scalar perturbations we achieve in this way 
azimuthal symmetry — the left hand side of the Boltz- 
mann equation and therefore also the brightness M de- 
pend only on fi — (k-n) and can be developed in Legendre 
polynomials. 

The left hand side of the Boltzmann equation for vector 
and tensor perturbations also determines the azimuthal 
dependence of M for vector and tensor perturbations, as 
we shall see in detail. 



Ms + ifikMs = 4z^A;($ - 
aaTTi, 



[D^;'^ -Ms^ 4zMl4 - \P2{f^)Q] , (65) 



Mf^ 



ifikMf^ 



-M 
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-(l-P2(/i))Q] 



(66) 



where 



Q = a. 
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The first term on the right hand side of Eq. ( |65| ) repre- 
sents the gravitational interaction (photons without col- 
lisions move along lightlike geodesies of the perturbed 
geometry), while the term in square brackets is the colli- 
sion integral for non-relativistic Compton scattering. 

Inserting expansion (Q) into Eqs. ( |65| ) and (|6^) using 
the standard recursion relations for Legendre polynomi- 
als, we obtain the following series of coupled equations: 



AS) 



ka 



, 



(67) 



■SS) 



(S) 



2a. 



^fc(*-$)-faaTne[^14-af^] , 



(68) 
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-[2a\ 



3(7. 



(69) 



^We could in principle add higher spin components to the 
distribution functions. But they are not seeded by gravity and 
since photons and neutrinos interact at high enough temper- 
atures, they are also absent in the initial conditions. 
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21+1 
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and 



for ^ > 3 



(70) 
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2£+l 
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( s) j_ 
+aaTne[-q^g + -Qi^m + -Se2)]- 



(71) 



For the neutrinos we obtain the same equations just 
without coUision integral 

oo 

AA,(n,k,t) = ^(-*)^(2^+l)z.f)(t,k)P,(A'), (72) 



and 



k 

21+1 
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(73) 



We are interested in the power spectrum of CMB 
anisotropies which is defined by 



ST 6T 
— (n) — (n) 



^^(2^+l)C,Pf(cos,?) 

(74) 
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Here (• • •) denotes the ensemble average over models. 
We assume that an 'ergodic hypothesis' is satisfied and 
we can interchange spatial and ensemble averages. The 
problem that actual observations can average at best over 
one horizon volume is known under the name 'cosmic 
variance'. It severely restricts the accuracy with which, 
for example, low multipoles of CMB anisotropies ob- 
served in our horizon volume can be predicted for a given 
model. 

Using the addition theorem of spherical harmonics, one 
obtains, with the Fourier transform conventions adopted 
here, (for details see Appendix B) 



riij 


= niTij - 






My 







riiTij - -Sij and 



We use coordinates for which k is parallel to the z-axis. 
Then 



(Si, £2,0) , uj^iu,.U2,0) 



and 



n = (i/l — /i^ cosif, — /i^ simp, fi) . 

With the ansatz 

7Wv(k,n,i) = Vl-M^[-^r^(k,Ai,t) cos(^- 
+7W^^^(k,^,t)sin^] , 



(77) 



the equations for A4i.2 decouple and the right hand side 
of Eq. (|7^) depends only on /x. Like for scalar perturba- 
tions, we expand M.1,2 in Legendre polynomials 



A^(^)(M,k,t) = ^(-i)^(2^+l)a(;;'(i,k)P,(/i), (78) 
£=0 



where e = 1, 2. 

Eq. ( |76| ) then leads to 
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With Eq. (jT^) , this can be expressed as the following set 
of coupled equations for the variables cr^V- 



^Z' + = aarn^lH'^ - a^,V] , (80) 



= ^ / k'dk{\<jf\to,k)\' 
on 



(75) 



where the superscript '^•^^ indicates that Eq. ( |75|) gives 
the contribution from scalar perturbations. 



2. vector perturbations 

Vector perturbations are very small on angular scales 
corresponding to ^ ^ 500, where Compton scattering and 
thus polarization become relevant. We therefore neglect 
polarization in this case. The Boltzmann equation for 
vector perturbations then reads 



Mv + ik • nMv = -4i(n • k)(n • S)- 
a(TTne[4(n • u)b) - Mv + ^ni-jMi-j] , 



(76) 



where 



a 



k 

e,l'- 3 

4 

^3 
and 



'e,0 



us- !^ (V) 1 (V) 

kY.^- aaTUel—crli — 



10 



10 



a 



(V) 



(V) 

2^-1-1 "" 

for£>2 



e,3 J ' 



KZi-(^+i)^Sli 



(81) 



(82) 



For neutrino perturbations we obtain the same equa- 
tions up to the collision term. We repeat them here for 
completeness. 



■,{V) 



21+1 



+ (83) 



As for scalar perturbations, the CMB anisotropy power 
spectrum is obtained by integration over fc-space. One 
finds (see Appendix B), 



7 



c. 



8tt 



(2£+l)2 



(84) 



Here the fact that there are two equal contributions from 
both polarization states, e — 1,2 (statistical isotropy) is 
taken care of. 



3. tensor perturbations 

For tensor perturbations, and a wave vector k pointing 
into the 3-direction, the only non vanishing components 
of the perturbed metric tensor are Hn = —H22 = 
and H12 — H21 = -ffx- Neglecting polarization, the 
Boltzmann equation for tensor perturbations is B 



Mt + ikfiMr — —An^n^ Hij 
-aarneiMr - ^n^jMij] . 
With the ansatz 

MT(k,n,t) = (1- fi'^)[M'^^\k,fi,t) cos 295 + 



7Wx^''(k, ^,t)sin2(p] , 



(85) 



(86) 



(T) 

the two modes ■A4\-i decouple completely and the right 
hand side of Eq. ( p5| ) depends only on /i. We can then 
expand the modes in terms of Legendre polynomials 

00 

M^P (m, k, t) = Y^{-if(n + iWS k)PHM), (87) 
where e = +, x. Eq. ( |8^ ) now becomes 

leading to the series of coupled equations for the coeffi- 
cients a 



(T) 



■JT) 

e,0 



I r (^) I (T) 
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10^ 



70 



■AT) 



21+1 



(T) 



for ^ > 1 . 



(89) 



(90) 



As before, the CMB anisotropy power spectrum is ob- 
tained by integration over fc-space (see Appendix B), 



(T) _ 1 (^ + 2)! 
87r {l--2)\ 



k^dk- 



{21+ If 



where 



a%l, 2{2£+l)aZ^ a^l, 
2i-l (2£-l)(2£ + 3) 2£ + 'S 



(91) 



(92) 



D. Eigenvector expansion of the source correlators 

In the previous subsections we have derived a closed 
system of linear differential equations with source terms. 
The source terms are linear combinations of the seed en- 
ergy momentum tensor which is determined by numerical 
simulations. A given realization of our model has random 
initial conditions; the seed energy momentum tensor is a 
random variable. In principle we could calculate the in- 
duced random variables L'g^'' (k, ig), K:(k,to), cr^*^(k,fo) 
etc for 100 to 1000 realizations of our model and de- 
termine the expectation values P{k), Pv{k) and Cn by 
averaging. This procedure has been adapted in Ref. p^ ] 
for a seed energy momentum tensor modeled by a few 
random parameters. 

In the case of a seed energy momentum tensor coming 
entirely from numerical simulations, this procedure is not 
feasible. The first and most important bottleneck is the 
dynamical range of the simulations which is about 40 
in our largest (400)'^ simulation, taking around 5 hours 
CPU time on a NEC SX-4 supercomputer. To determine 
the Cf's for 2 < ^ < 1000 we need a dynamical range of 
about lO'OOO in fc-space (this means fcmax/^min ^ lO'OOO, 
where fcmax and fcmin are the maximum and minimum 
wave numbers which contribute to the C^'s within our 
accuracy ( ~ 10%). 

With brute force, this problem is thus not tractable 
with present or near future computing capabilities. But 
there are a series of theoretical observations which reduce 
the problem to a feasible one: 

For each wave vector k given, we have to solve a system 
of linear perturbation equations with random sources. 



VX 



(93) 



Here I? is a time dependent linear differential operator, 
X is the vector of our matter perturbation variables spec- 
ified in the previous subsections (photons, CDM, baryons 
and neutrini; total length up to 2000), and S is the ran- 
dom source term, consisting of linear combinations of the 
seed energy momentum tensor. 

For given initial conditions, this equation can be solved 
by means of a Green's function (kernel), Q{t,t'), in the 
form 



X, (to, k) = ^ dtg^i {to, t, k)5i (i, k) . 



(94) 



We want to compute power spectra or, more generally, 
quadratic expectation values of the form 

(x,(to,k)xr(to,k)) , 

which, according to Eq. (|9j) are given by 

(x,(to,k)xr(to,k)) = 

tin Jtin 



(95) 
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The only information about the source random variable 
which we really need in order to compute power spectra 
are therefore the unequal time two point correlators 

{SUt,k)S:{t'M . (96) 

This nearly trivial fact has been exploited by many work- 
ers in the field, for the first time probably in Ref. p^ ] 
where the decoherence of models with seeds has been 
discovered, and later in Refs. P jl9|j20| , [T^ and others. 

To solve the enormous problem of dynamical range, we 
make use of 'scaling', statistical isotropy and causality. 

We call seeds 'scaling' if their correlation functions 
Cfj,i,px defined by 

Q^,,{k,t) = M^e^,{k,t) , (97) 
C^,pA(k,i,t') = (0f,4k,t)0;^{k,t')) (98) 

are scale free; i.e. the only dimensional parameters 
in CpjypA are the variables t, t' and k themselves. Up 
to a certain number of dimensionless functions F„ of 
z = kVtF and r — t/t' , the correlation functions are then 
determined by the requirement of statistical isotropy, 
symmetries and by their dimension. Causality requires 
the functions F„ to be analytic in z^. A more detailed 
investigation of these arguments and their consequences 
is presented in Ref. There we show that statistical 
isotropy and energy momentum conservation reduce the 
correlators ( p8[ ) to five such functions Fi to F^. 

In cosmic string simulations, energy and momentum 
are not conserved. Strings loose their energy by radiation 
of gravitational waves and/or massive particles. In this 
case 14 functions of and r are needed to describe the 
unequal time correlators ||2l|] . 

Since analytic functions generically are constant for 
small arguments <C 1, i^„(0,r) actually determines F„ 

for all values of k with z = k^/tF ~ 0.5. Furthermore, 
the correlation functions decay inside the horizon and we 
can safely set them to zero for z ~ 40 where they have 
decayed by about two orders of magnitude (see Figs. 1 
to 11). Making use of these generic properties of the 
correlators, we have reduced the dynamical range needed 
for our computation to about 40, which can be attained 
with the (256)'^ to (512)^ simulations feasible on present 
supercomputers. 

For the scalar part we need the correlators 

(<i>,(k,i)cl>:(k,i')) = -^j^Cn(z,r) , (99) 

($,(k,t)vl':(k,0) = -^^Ci2(z,r) , (100) 

(*,(k,t)vl/:(k,0) = -^C22(z,r) , (101) 

as well as C2i(z,r) — Ci2{z, ^/r). The functions Cij are 
analytic in z^. The pre- factor \/ (k^^fW) comes from the 
fact that the correlation functions {fpfp), k'^{J-„f*) and 



ifvfv) have to be analytic and from dimensional consid- 
erations (see Ref. |^). 

The functions Cij are shown in Figs. 1 to 3. Panels (a) 
are obtained from numerical simulations. Panels (b) rep- 
resent the same correlators for the large- TV limit of global 
0(A^)-models (see |2|,|§). 




FIG. 1. The two point 

correlation function Cii(z,r) = fc'' Vtt^{'l>s (k, i)<l>* (k, i')) is 
shown. Panel (a) represents the result from numerical sim- 
ulations of the texture model; panel (b) shows the large- 
limit. For fixed r the correlator is constant for z < 1 and 
then decays. Note also the symmetry under r — + 1/r. 

In Fig. 4 we show Cij{z,r = 1), and in Fig. 5 the 
'constant' of the Taylor expansion for dj is given as a 
function of r, i.e., C.y(0,r). 

Vector perturbations are induced by S which is 
seeded by w^'"\ Transversality and dimensional argu- 
ments require the correlation function to be of the form 

{w'i^\k,t)w'^^^*ik,t')) = VtPik^d^, - k,k,)W{z,r) . 

(102) 

Again, as a consequence of causality, the function W is 
analytic in z^ (see [Q). The function W{z,r) is plotted 
in Fig. 6. In Figs. 7 and 8 we graph W{z, 1) and W{0, r). 
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FIG. 2. 



0.01 



but for 



Symmetry, transversality and tracelessness, together 
with statistical isotropy require the tensor correlator to 
be of the form (see Q]) 



:T{z, r)[SiiSjjn + SimSjl ~ (5.y (J/m + k [Sijkihm + 



^imkikj Silkjkjji Simkikj Sjikikjji Sjtnkiki) + 
k kikjkikfYi^ ■ (103) 

The functions T{z,r) as well as T(z, 1) and T{0,r) are 
shown in Figs. 9 to 11. 

Clearly, all correlations between scalar and vector, 
scalar and tensor as well as vector and tensor pertur- 
bations have to vanish. 

The scalar source correlation matrix C and the func- 
tions W and T can be considered as kernels of positive 
hermitian operators in the variables x = kt = zr^l"^ and 
x' = kt' = z/r^/^, which can be diagonalized. 



a,{x,x')=Y.^nMn\^)vS^*i^'), 

71 

Wix,x') = Y,X(r\(y{x)v(J^*{x'), 



nx,x') = Y.^(j\(p{x)^ 



XT) 



(104) 
(105) 
(106) 




0.01 



FIG. 3. The unequal time correla- 

tor, \Cx2{,z,r)\ = A;''ytt^|(<I>4k,t)**(k,t'))j is shown. Note 
that the r ^ 1/r symmetry is lost in this case. 

where the series ^^'^ (^"^0 

mal series of eigenvectors (ordered according to the am- 
plitude of the corresponding eigenvalue) of the operators 
C, W and T respectively for a given weight function w. 
We then havePI 



a,ix,x')vfj{x')w{x')dx' = X^hl^ix) , 
Wix,x')vl^Hx')wix')dx' = X'^/^vl^Hx) 
T{x,x')v^P{x')w{x')dx' - )«f )(x) 



(107) 
(108) 
(109) 



The eigenvectors and eigenvalues depend on the weight 
function w which can be chos en t o op timize the speed of 
convergence of the sums (104) to (106). In our models we 
found that scalar perturbations typically need 20 eigen- 
vectors whereas vector and tensor perturbations need five 



^Here the assumption that the operators C, W and T are 
trace-class enters. This hypothesis is verified numerically by 



the fast convergence of the sums (104) to (106) 



10 



0.1 



0.01 0.1 1 
z 



10 



0.1 



_l I 



10 



100 



I I 1 1 Mill — I I I iiiii| — I I 1 1 iiiij — ri 




FIG. 4. The correlators Cij{z,l) are shown. The solid, 
dashed and dotted hnes represent C22 , Cn and iCi2| respec- 
tively. Panel (a) is obtained from numerical simulations of the 
texture model and panel (b) shows the large- A'' limit. A strik- 
ing difference is that the large- A'^ value for IC12I is relatively 
well approximated by the perfectly coherent result y^\CiiC22\ 
while the texture curve for IC12I Hes nearly a factor 10 lower. 



to ten eigenvectors for an accuracy of a few percent (see 
Fig. [1]). 

Inserting Eqs. ( |lO^ to ( [l06| ) in Eq. (H), leads to 



(X,(k,io)^;(k,to)) -^A„xr(fcio)xj")*(Ho) , 



(110) 



where X^^\to) is the sohition of Eq. (|93|) with determin- 



istic source term v, 



(") 



xj")(io,k) = ^ " dtg{to,t,k)jivl"\x,k) . (Ill) 
For the CMB anisotropy spectrum this gives 

ns riv fiT 

n n n 

(112) 

cl'"-* is the CMB anisotropy induced by the determinis- 
tic source w„, and n, is the number of eigenvalues which 
have to be considered to achieve good accuracy. 



100 



FIG. 5. The correlators C, 




are shown in the same 



line styles as in Fig. ^, but for z — as function of r = t' /t. 
The stronger decoherence of the texture model is even more 
evident here. 



Instead of averaging over random solutions of Eq. (p4|), 
we can thus integrate Eq. (^4|) with the deterministic 
source term w'"^ and sum up the resulting power spectra. 
The computational requirement for the determination of 
the power spectra of one seed model with given source 
term is thus on the order of ns inflationary models. This 
eigenvector method has first been applied in Ref. |^ . 

A source is called totally coherent |^,^ if the unequal 
time correlation functions can be factorized. This means 
that only one eigenvector is relevant. A simple totally 
coherent approximation, which however misses some im- 
portant characteristics of defect models, can be obtained 
by replacing the correlation matrix by the square root of 
the product of equal time correlators, 



{ms;it')) ±^{\sdm{\sAt'W) • (113) 

This approximation is exact if the source evolution is lin- 
ear. Then the different k modes do not mix and the value 
of the source term at fixed k at a later time is given by its 
value at initial time multiplied by some transfer function, 
S(k, t) — 5(k, ti„)r(k, t, tin). In this situation, (113) be- 
comes an equality and the model is perfectly coherent. 
Decoherence is due to the non-linearity of the source evo- 
lution which induces a 'sweeping' of power from one scale 
into another. Different wave numbers k do not evolve in- 
dependently. 

It is interesting to note that the perfectly coherent ap- 
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FIG. 6. The vector correlator W{z,r) is shown. The 
texture simulations, panel (a), and the large- A'^ limit, panel 
(b), give very similar results. 



FIG. 7. The vector correlator W{z, 1) is plotted. The solid 
line represents the texture simulations and the dashed line is 
the large- result. Up to a slight difference in amplitude, the 
two results are very similar. 
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FIG. 8. The vector correlator W(0,r) is shown. The solid 
line represents the texture simulations and the dashed line is 
the large- A'' result. Also here, the two results are very similar. 
The 'wings' visible in the texture curve are probably not due 
to a resolution problem but the beginning of oscillations. 



proximation, (113), leaves open a choice of sign which has 
to be positive if i = j, but which is undetermined oth- 
erwise. According to Schwarz inequality the correlator 
(5,(<)5*(f )) is bounded by 



y^(|5,(i)P)(|5,(f)P)<(5,:(t)5;(0) < 
(|5.(t)P)(|5,(t')P). (114) 



Therefore, for scales/variables for which the Greens func- 
tion is not oscillating (e.g. Sachs Wolfe scales) the full 
result always lies between the 'anti-coherent' (minus sign) 
and the coherent result. We have verified this behavior 
numerically. 

The first evidence that Doppler peaks are suppressed 
in defect models has been obtained in the perfectly co- 
herent approximation in Ref . . In Fig. |l3| we show the 
contributions to the C^'s from more and more eigenvec- 
tors. A perfectly coherent model has only one non-zero 
eigenvalue. 

A comparison of the full result with the totally coher- 



sees that decoherence does smear out the oscillations 
present in the fully coherent approximation, and does 
somewhat damp the amplitude. Decoherence thus pre- 
vents the appearance of a series of acoustic peaks. The 
absence of power on this angular scale, however, is not 
a consequence of decoherence but is mainly due to the 
anisotropic stresses of the source which lead to pertur- 
bations in the geometry inducing large scale C^'s (Sachs 
Wolfe), but not to density fluctuations. Large anisotropic 
stresses are also at the origin of vector and tensor fluctu- 
ations. Our results are in agreement with Refs. and 
Q but we disagree with Ref. , which has found acous- 
tic peaks with an amplitude of about six in the coherent 
approximation. 

In the real universe, perfect scaling of the seed corre- 
lation functions is broken by the radiation-matter tran- 
sition, which takes place at the time of equal matter and 



radiation, t 



7 ''eg 



ent approximation is presented in Fig. 14, There one 



20/i"^rj„ ' Mpc. The time ieg IS an 
additional scale which enters the problem and influences 
the seed correlators. Only in a purely radiation or mat- 
ter dominated universe are the correlators strictly scale 
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FIG. 9. As Fig. g, but for the tensor source function 
T{z,r). 

invariant. This means actually that the k dependence of 
the correlators C, W and T cannot really be cast into 
a dependence on x and x' , but that these functions de- 
pend on t, t' and fc in a more complicated way. We have 
to calculate and diagonalize the seed correlators for each 
wave number k separately and the huge gain of dynami- 
cal range is lost as soon as scaling is lost. 

In the actual case at hand, however, the deviation from 
scaling is weak, and most of the scales of interest to us 
enter the horizon only in the matter dominated regime. 
The behavior of the correlators in the radiation domi- 
nated era is of minor importance. To solve the problem, 
we calculate the correlator eigenvalues and eigenfunctions 
twice, in a pure radiation and in a pure matter universe 
and we interpolate the source term from the radiation to 
the matter epoch. Denoting by Am, Wm and Xr,Vr a given 
pair of eigenvalue and eigenvector in a matter and radia- 
tion universe respectively, we choose as our deterministic 
source function 

v{t) = y{t)y/Kvr{kt) + (1 - y{t))y/X^Vm{kt) (115) 
with, e.^.. 
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FIG. 10. As Fig. ^, but for the tensor source function 
T{z,l). 
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r 

FIG. 11. As Fig. ^, but for the tensor source function 

r(o,r). 

y{t) = 7^ or y{t) = eM-t/te,) , (116) 

or some other suitable interpolation function. In Fig. |l^ 
we show the results for scalar, vector and tensor pertur- 
bations respectively using purely radiation dominated era 
and from interpolated source terms. 

Clearly the effect of the radiation dominated early 
state of the universe is relatively unimportant for the 
scales considered here. The difference between the pure 
matter era result and the interpolation is barely visible 
and thus not shown on the plot. This seems to be quite 
different for cosmic strings where the fluctuations in the 
radiation era are about twice as large as those in the mat- 
ter era The radiation dominated era has very little 
effect on the key results which we are reporting here; 
namely the absence of acoustic peaks and the missing 
power on very large scales. 

In models with cosmological constant, there is actually 
a second break of scale invariance at the matter-A tran- 
sition. There we proceed in the same way as outlined 
above. Since defects cease to scale and disappear rapidly 
in an exponentially expanding universe, the eigenvalues 
for the A dominated universe all vanish. 
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FIG. 12. The sum of the first few eigenfunctions of T{x, x) 
is shown for two different weight functions, (a) logarithmic, 
w — 1/x and (b) hnear, lo = 1. The first (long dashed), first 
and second (short dashed), first ten (dotted) and first thirty 
(solid) eigenfunctions are summed up. The open circles repre- 
sent the full correlation function. Clearly, the eigenfunctions 
obtained by linear weighting converge much faster. Ifere we 
only show the equal time diagonal of the correlation matrix, 
but the same behavior is also found in the Ce power spectrum 
which is sensitive to the full correlation matrix. 

E. Initial conditions and numerical implementation 

We numerically integrate our system of equations 
from redshift z — 10^ up to the present with the goal 
to have one percent accuracy up to £ '-^ 1000, for a 
given source term. We use the integration method de- 
scribed in Refs. jl^ and [|^. We sample the interval 
—5 < logj^Q fc/i~^Mpc < —0.75 with minimum step size 
Alogj^Q = 0.04, for the scalar case and use a smooth- 
ing algorithm to suppress high frequency sampling noise. 
In order to save computing time, we start the inte- 
gration of the (T^(fc)'s with 10 harmonics, adding new 
harmonics in the course of the integration. We find 
that typically ^ 40 harmonics are sufficient for small 
k values (logj^g A:/i~^Mpc ^ —3), while for higher k 
(logj^g A:/i~"'^Mpc ^ —1), up to ^ 1500 harmonics for 
the scalar case, ~ 200 for the tensor case are needed to 
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FIG. 13. The scalar, vector and tensor contributions for 
the texture model of structure formation are shown. The 
dashed lines show the contributions from single eigenfunctions 
while the solid line represents the sum. Note that the single 
contributions to the scalar and tensor spectrum do show os- 
cillations which are however washed out in the sum. (Vector 
perturbations do not obey a wave equation and thus do not 
show oscillations.) 



achieve the desired accuracy. Including more than 40 har- 
monics for neutrinos corrects our results by less than one 
percent. We obtain <i> algebraically using Eq. (p^. With 
this choice of variables we avoid the numerical difficulties 
present in conformal gauge pst , where $ is determined 
by numerical integration. 

The abundance of free electrons, Ue, is calculated fol- 
lowing a standard recombination scheme j29j for H and 
^He, for a helium abundance by mass of 23%. At high 
redshift z > 10^, the Thomson opacity is very large, 
and photons and baryons are tightly coupled. Due to 
the large Thomson drag term, Eqs.(^) and ( p8| ) become 
stiff and difficult to solve numerically. Therefore, in this 
limit we follow the method of Ref . 1 30 1 , which is accurate 
to second order in (crT"-e)~"'^ (see also [^). Assuming 
a standard inflationary model, we obtain a single scalar 
power spectrum in few minutes (^ 30 seconds for the 
tensor case) on a PC class workstation, which differs by 
less than one percent from the Cf's computed with other 
codes [0, m. 

Summing the scalar C^'s from the largest 15 eigenvec- 
tors (5 in the tensor case, 10 for vector perturbations) 
typically reproduces the total sum to better than 5% (see 
Fig. IH). 
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FIG. 14. The Ci power spectrum for the texture scenario 
is shown in the perfectly coherent approximation (top panel) 
and in the full eigenfunction expansion. Even in the coher- 
ent approximation, the acoustic peaks are not higher than 
the Sachs Wolfe plateau. Decoherence just washes out the 
structure but does not significantly damp the peaks. 



- Including radiation to matter transition 

- Radiation dominated 




I 

FIG. 15. The scalar, vector, tensor and total Ce power 
spectrum is shown from pure radiation sources and for an in- 
terpolated source. While the vector perturbations are some- 
what higher in the radiation era, scalar and tensor pertur- 
bations are higher in the matter era and the sum is nearly 
unchanged. 



III. THE NUMERICAL SIMULATIONS 

As in previous work [Q, we consider a spontaneously 
broken scalar field with O(A^) symmetry. We use the 
CT- model approximation, i.e., the equation of motion 



□/3-(/3-n/3)/3 = 0, 



(117) 



where j3 is the rescaled field /3 — (p/rj. 

We do not solve the equation of motion directly, but 
use a discretized version of the action Ip2i : 



A 



d^(3 • + - - 1) 



, (118) 



where A is a Lagrange multiplier which fixes the field 
to the vacuum manifold (this corresponds to an infinite 
Higgs mass) . Tests have shown that this formalism agrees 
well with the complementary approach of using the equa- 
tion of motion of a scalar field with Mexican hat potential 
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FIG. 16. 

C22(Z,1) = 
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kt 

The ETCs Cii{z,l) = (l$P)(fct) (panel a) and 
1*1^) (fct) (panel b) are shown for different times. 
In grid units the times are t — 4 (dashed), t = 8 (dotted), 
t — 12 (long dashed), t = 16, 20 (dash dotted, long dash 
dotted) and t = 24 (solid). Clearly C22 scales much sooner 
than Cii. To safely arrive in the scaling regime one has to 
wait until t ~ 16 and Cij{kt = 0) is best determined at t > 20 
but kt < 1. 



and setting the inverse mass of the particle to the smallest 
scale that can be resolved in the simulation (typically of 
the order of 10~^^ GeV), but tends to give better energy 
momentum conservation. 

As we cannot trace the field evolution from the unbro- 
ken phase through the phase transition due to the limited 
dynamical range, we choose initially a random field at a 
comoving time t — 2Ax. Different grid points arc uncor- 
related at all earlier times |Q . 

The use of finite differences in the discretized action 
as well as in the calculation of the energy momentum 
tensor introduce immediately strong correlations between 
neighboring grid points. This problem manifests itself in 
an initial phase of non-scaling behaviour, the length of 
which varies between lOAa; and 20Ax, depending on the 
variable considered. It is very important to use results 
from the scaling regime only (cf. Fig. [l^ ). 

In order to reduce the time necessary to reach scaling 
and to improve the overall accuracy, we try to choose the 
finite differences in an optimal way. Our current code 
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calculates all values in the center of each cubic cell defined 
by the lattice. The additional smoothing introduced by 
this improves energy-momentum conservation by several 
percent.^ 

To calculate the unequal time correlator (UTC), the 
value of the observable under consideration is saved once 
scaling is reached at time t^. (we checked this by using dif- 
ferent correlation times) and then correlated at all follow- 
ing time steps. While there is some danger of contaminat- 
ing the equal time correlator (ETC), which contributes 
most strongly to the C^'s, with non-scaling sources, this 
method ensures that the constant for kt ^ is deter- 
mined with maximal precision for the ETCs. This is 
very important as the constants Cij(0,l), 1^(0,1) and 
T(0, 1) fix the relative size of scalar, vector and tensor 
contributions of the Sachs- Wolfe part and severely influ- 
ence the resulting Q's. In contrast, the CMB spectrum 
seems quite stable under small variations of the shape of 
the UTCs. 

The resulting UTCs are obtained numerically as func- 
tions of the variables k, t and with t > tc and tc fixed. 
They are then linearly interpolated to the required range. 
We construct a hermitian 100 x 100 matrix in kt and kt', 
with the values of kt chosen on a linear scale to maximize 
the information content, < fci < Xmax- The choice 
of a linear scale ensures good convergence of the sum 
of the eigenvectors after diagonalization (see Fig. 12), 



but still retains enough data points in the critical region, 
0{x) = 1, where the correlators start to decay. In prac- 
tice we choose as the endpoint Xmax of the range sampled 
by the simulation the value at which the correlator decays 
by about two orders of magnitude, typically Xmax ~ 40. 
The eigenvectors that are fed into the Boltzmann code 
are then interpolated using cubic splines with the condi- 
tion Vn{kt) for kt ^ a^max- 

We use several methods to test the accuracy of the 
simulation: energy momentum conservation of the de- 
fects code is found to be better than 10% on all scales 
larger than about 4 grid units, as is seen in Fig. A 
comparison with the exact spherically symmetric solution 
in non-expanding space [ p2[ shows very good agreement. 

The resulting CMB spectrum on Sachs Wolfe scales is 
consistent with the line of sight integration of Ref. [Q. 
Furthermore, the overall shape and amplitude of the un- 
equal time correlators are quite similar to those found in 
the analytic large- iV approximation (see Figs. 1 

to 11). The main difference of the large- iV approxima- 



^ Julian Borrill suggested to introduce "spherical derivatives" 
that take into account the fact that the vacuum manifold is a 
N-sphere and therefore curved, and that this curvature should 
be important at least in the initial stages of the simulation 
and for unwinding events . So far we haven't investigated 
this idea sufficiently to include it into our production code. 
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FIG. 17. Energy momentum conservation of our numeri- 
cal simulations is shown. The lines represent the sum of the 
terms which has to vanish if energy (solid) respectively mo- 
mentum (dashed) is conserved, divided by the sum of the ab- 
solute value of these terms. The abscissa indicates the wave- 
length of the perturbation as fraction of the size of the entire 
grid. 



tion is that there the field evolution, Eq. (117), is ap- 
proximated by a linear equation. The non-linearities in 
the large- seeds which are due solely to the energy mo- 
mentum tensor being quadratic in the fields, are much 
weaker than in the texture model where the field evolu- 
tion itself is non-linear. Therefore, decoherence which is 
a purely non-linear effect, is expected to be much weaker 
in the large- iV limit. This is actually the main difference 
between the two models as can be seen in Fig. nsl 



IV. RESULTS AND COMPARISON WITH DATA 

A. CMB anisotropies 

The Cf's for the 'standard' global texture model are 
shown in Fig. |lj, bottom panel. 
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TABLE I. The value of the normalization constant e and 
the fluctuation amplitude erg are given for the different mod- 
els considered. The error in e comes from a best fit normal- 
ization to the full CMB data set. Cosmological parameters 
which are not indicated are identical in all models or given 
by fio = f^cdm + SI A -I- Sib = 1. We consider only spatially fiat 
models with fi;, = 0.05 and a helium fraction of 23%. The 
parameter choice indicated in the top line is referred to as 
standard texture model in the text. 
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FIG. 18. Top panel: the — model. Bottom panel: 
The Ce power spectrum is shown for the large- A'' limit (bold 
line) and for the texture model. The main difference is clearly 
that the large- A'^ curve shows some acoustic oscillations which 
are nearly entirely washed out in the texture case. 
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FIG. 19. The Ci power spectrum is shown for different 
values of cosmological parameters. In the top panel we choose 
JIa ~ 0, QcDM ~ 0.95, Qb = 0.05 and vary h. In the bottom 
panel we fix ft = 0.5, fii, — 0.05 and vary Qa- We only 
consider spatially flat universes, Oo = 1. 



Vector and tensor modes are found to be of the same 
order as the scalar component at COBE-scales. For the 



'standard' texture model we obtain c[q^ : C'^n ^ : C- 



^10 



'10 



0.9 : 1.0 : 0.3, in good agreement with the predictions of 
Refs. ||Jl^,|l^ and Due to tensor and vector con- 
tributions, even assuming perfect coherence (see Fig. |l^, 
top panel), the total power spectrum does not increase 
from large to small scales. Decoherence leads to smooth- 
ing of oscillations in the power spectrum at small scales 
and the final power spectrum has a smooth shape with 
a broad, low isocurvature 'hump' at £ ~ 100 and a small 
residual of the first acoustic peak at £ ~ 350. There is no 
structure of peaks at small scales. The power spectrum 
is well fitted by the following fourth-order polynomial in 
X = log£: 



HOC 



1.5 - 2.6x + 3.3x^ - lAx-^ + 0.17x'^ 



10 



(119) 



The effect of decoherence is less important for the 
large- model, where oscillations and peaks are still vis- 
ible (see Fig ^ bottom panel). This is due to the 
fact that the non-linearity of the large- limit is only 
in the quadratic energy momentum tensor. The scalar 
field evolution is linear in this limit |^^, in contrast to 
the N = 4 texture model. Since decoherence is inher- 
ently due to non-linearities, we expect it to be stronger 
for lower values of A^. COBE normalization leads to 
e = (0.92 ± 0.1)10-5 for the large- A^ limit. 

In Fig. ^ we plot the global texture Ci power spectrum 
for different choices of cosmological parameters. The 
variation of parameters leads to similar effects like in the 
inflationary case, but with smaller amplitude. At small 
scales {£ > 200), the Q's tend to decrease with increas- 
ing Hq and they increase when a cosmological constant 
= 1 — fim is introduced. Nonetheless, the amplitude 
of the anisotropy power spectrum at high ^s remains in 
all cases on the same level like the one at low is, without 
showing the substantial peak found in inflationary mod- 
els. The absence of acoustic peaks is a stable prediction 
of global 0{N) models. The models are normalized to 
the full CMB data set, which leads to slightly larger val- 
ues of the normalization parameter e = AnGr]^ than pure 
COBE normalization. In Table | we give the cosmologi- 
cal parameters and the value of e for the models shown 
in Fig. |l|. 

In order to compare our results with current exper- 
imental data, we have selected a set of 31 different 
anisotropy detections obtained by different experiments, 
or by the same experiment with different window func- 
tions and/or at different frequencies. Theoretical predic- 
tions and data of CMB anisotropics are usually compared 
by plotting the theoretical Cg curve along with the CMB 
measurements converted to band power estimates. We do 
this in the top panel of Fig. |2^. The data points show an 
increase in the anisotropics from large to smaller scales. 
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FIG. 20. The Ct spectrum obtained in the standard tex- 
ture model is compared with data. In the top panel exper- 
imental results and the theoretical curve are shown as func- 
tions of i. In the two lower panels we indicate the value of 
each of the 31 experimental data points with 1-a error bars 
and the corresponding theoretical value with its uncertainty. 
The experiments corresponding to a given number are given 
in Table In the middle panel the 8 COBE data points are 
shown. In the bottom panel other experiments are presented. 



in contrast to the theoretical predictions of the model. 
This fashion of presenting the data is surely correct, but 
lacks informations about the uncertainties in the theoret- 
ical model. Therefore we also compare the detected mean 
square anisotropy, A*^-^^^'^ and the experimental 1-a er- 
ror, I](^^p)2^ directly with the corresponding theoretical 
mean square anisotropy, given by 



A^^'*) = ^Y,{2i + l)CiWi 



(120) 



where the window function Wi contains all experimental 
details (chop, modulation, beam, etc.) of the experiment. 

The theoretical error in principle depends on the statis- 
tics of the perturbations. If the distribution is Gaussian, 
one can associate a sample/cosmic variance 



1 1 

J8^ 



(121) 



where / represents the fraction of the sky sampled by a 
given experiment. 



Deviation from Gaussianity leads to an enhancement 
of this variance, which can be as large as a factor of 
7 (see iQ). Even if the perturbations are close to 
Gaussian (which has been found by simulations on large 
scales |7|j37|]), the Q's, which arc the squares of Gaus- 
sian variables, are non-Gaussian. This effect is, however 
only relevant for relatively low ^s. Keeping this caveat 
in mind, and missing a more precise alternative, we nev- 
ertheless indi cate t he minimal, Gaussian error calculated 
according to ( 121 ). We add a 30% error from the CMB 
normalization. The numerical seeds are assumed to be 
about 10% accurate. 

In Table || the detected mean square anisotropy, 
A(-^^p), with the experimental l-tr error are listed for 
each experiment of our data set. The corresponding sky 
coverage is also indicated. In Fig. |2^ we plot these data 
points, together with the theoretical predictions for a tex- 
ture model with h = 0.5 and I^a = 0. 

We find that, apart from the COBE quadrupole, only 
the Saskatoon experiment disagrees significantly, more 
than 1(7, with our model. But also this disagreement is 
below 3(7 and thus not sufficient to rule out the model. 
In the last column of Table ^ we indicate 



(A} 



(Th) 



(Exp) ^2 



)7(E 



{Th)2 
j 



for the j-th experiment, where the theoretical model 
is the standard texture model with = and h = 
0.5. The major discrepancy between data and the- 
ory comes from the COBE quadrupole. Leaving away 
the quadrupole, which can be contaminated and leads 
to a similar also for inflationary models, the data 
agrees quite well with the model, with the exception 
of three Saskatoon data points. Making a rough chi- 
square analysis, we obtain (excluding the quadrupole) a 
value — J2j x] 30 for a total of 30 data points 
and one constraint. An absolutely reasonable value, but 
one should take into account that the experimental data 
points which we are considering are not fully indepen- 
dent. The regions of sky sampled by the Saskatoon 
and MSAM or COBE and Tenerife, for instance, over- 
lap. Nonetheless, even reducing the degrees of freedom 
of our analysis to iV = 25, our is still in the range 
(TV - 1) ± v^(7V - 1) ~ 24 ± 7 and hence still compatible 
with the data. 

This shows that even assuming Gaussian statistics, the 
models are not convincingly ruled out from present CMB 
data. There is however one caveat in this analysis: A 
chi-square test is not sensitive to the sign of the discrep- 
ancy between theory and experiment. For our models 
the theoretical curve is systematically lower than the ex- 
periments. For example, whenever the discrepancy be- 
tween theory and data is larger than 0.5tT, which hap- 
pens with nearly half of the data points (13), in all cases 
except for the COBE quadrupole, the theoretical value is 
smaller than the data. If smaller and larger are equally 
likely, the probability to have 12 or more equal signs is 
2(13 + l)/2i3 ~ 3.4 X 10-^. This indicates that either 
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the model is too low or that the data points are system- 
atically too high. The number 0.003 can however not 
be taken seriously, because we can easily change it by 
increasing our normalization on a moderate cost of x^. 



B. Matter distribution 

In Table || we show the expected variance of the total 
mass fluctuation an in a ball of radius R — S/i^^Mpc, 
for different choices of cosmological parameters. We 
find (78 ~ (0.44 ± 0.07) ft. (the error coming from the 
CMB normalization) for a flat model without cosmolog- 
ical constant, in agreement with the results of Ref. 
From the observed cluster abundance, one infers (jg ~ 
(0.50 ± 0.04)^2-0-5 g ^ 0.59t[! j^ |l[. These 

results, which are obtained with the Press-Schechter for- 
mula, assume Gaussian statistics. We thus have to take 
them with a grain of salt, since we do not know how non- 
Gaussian fluctuations on cluster scales arc in the texture 
model. According to Ref. |Q, the Hubble constant lies 
in the interval h ~ 0.73 ± 0.06 ± 0.08. Hence, in a flat 
CDM cosmology, taking into account the uncertainty of 
the Hubble constant, the texture scenario predicts a rea- 
sonably consistent value of (Ts- 

As already noticed in Refs. |l^ and [||, unbiased global 
texture models are unable to reproduce the power of 
galaxy clustering at very large scales, ~ 2Qh^^ Mpc. In 
order to quantify this discrepancy we compare our pre- 
diction of the linear matter power spectrum with the re- 
sults from a number of infrared ( | 



^■^-0.15 



|54| ) and optically- 
selected ( [^) galaxy redshift surveys, and with the 
real-space power spectrum inferred from the APM pho- 
tometric sample ( Q) (see Fig. |2^). Here, cosmological 
parameters have important effects on the shape and am- 
plitude of the matter power spectrum. Increasing the 
Hubble constant shifts the peak of the power spectrum 
to smaller scales (in units of ft/Mpc), while the inclusion 
of a cosmological constant enhances large scale power. 

We consider a set of models in f^A - h space, with linear 
bias ||53| as additional parameter. In Table 111 we report 



for each survey and for each model the best value of the 
bias parameter obtained by x^-niinimization. We also 
indicate the value of (not divided by the number of 
data points). The data points and the theoretical predic- 
tions are plotted in Fig. 21. Our bias parameter strongly 
depends on the data considered. This is not surprising, 
since also the catalogs are biased relative to each other. 

Models without cosmological constant and with h ^ 
0.8 only require a relatively modest bias 6 ~ 1.3 — 3. 
But for these models the shape of the power spectrum 
is wrong as can be seen from the value of which is 
much too large. The bias factor is in agreement with our 
prediction for ug. For example, our best fit for the IRAS 
data, for /i - 0.8 is & ~ 1.3. With cr|^^^ = (0.69 ±0.05), 
this gives erg ~ 0.48 ± 0.04, compatible with the direct 
computation 



Whether IRAS galaxies are biased is still under debate. 
Published values for the /3 parameter, defined as /3 

for IRAS galaxies, range between /?/ 
|9| and (3i = 0.5 ±0.1 g]. Biasing of IRAS galaxies 
is also suggested by measurements of bias in the optical 
band. For example, Ref. ||6l| finds j3o — 0.40 ± 0.12, in 
marginal agreement with [ |62[ , which obtains /3o ~ 0.35 ± 
0.1. A bias for IRAS galaxies is not only possible but 
even preferred in flat global texture models. 

But also with bias, our models are in significant con- 
tradiction with the shape of the power spectrum at large 
scales. As the values of in Table HI and Fig. ^ clearly 
indicate, the models are inconsistent with the shape of 
the IRAS power spectrum, and they can be rejected with 
a high confidence level. The APM data which has the 
smallest error bars is the most stringent evidence against 
texture models. Nonetheless, these data points are not 
measured in redshift space but they come from a de- 
projection of a 2 — D catalog into 3 — D space. This 
might introduce systematic errors and thus the errors of 
APM may be underestimated. 

Models with a cosmological constant agree much better 
with the shape of the observed power spectra, the value 
of being low for all except the APM data. But the 
values of the bias factors are extremely high for these 
models. For example, IRAS galaxies should have a bias 
h ^ 3 — 6, resulting in erg < 0.25, and in a /?/ < 0.2 
which is too small, even allowing for big variances due to 
non-Gaussian statistics. 

The power spectra for the large- iV limit and for the co- 
herent approximation are typically a factor 2 to 3 higher 
(see Fig. 22), and the biasing problem is alleviated for 
these cases. For JIa = we find cg = 0.57ft. for the large- 
N limit and erg — 0.94ft. for the coherent approximation. 
This is no surprise since only one source function, ^^s, 
the analog of the Newtonian potential, seeds dark matter 
fluctuations and thus the coherence always enhance s the 
unequal time correlator. The second inequality in (114) 
applies. The dark matter Greens function is not oscil- 
lating, so this enhancement translates directly into the 
power spectrum. 

Models which arc anti-coherent in the sense defined in 
Section IID reduce power on Sachs- Wolfe scales and en- 
hance the power in the dark matter. Anti-coherent scal- 
ing seeds arc thus the most promising candidates which 
may cure some of the problems of global 0{N) models. 

The simple analysis carried out here does not take into 
account the effects of non-linearities and redshift distor- 
tions. Redshift distortions in the texture case should 
be less important than in the inflationary case since the 
peculiar velocities are rather low (see next paragraph). 
Non-linearities typically set in at A; > 0.5ftMpc-^ and 
should not have a big effect on our main conclusions 
which come from much larger scales. Inclusion of these 
corrections will result in more small-scale power and in 
a broadening of the spectra, which even enhances the 
conflict between models and data. Furthermore, varia- 
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FIG. 21. The dark matter power spectrum for the texture 
model (solid line) is compared with the coherent approxima- 
tion (short dashed) and the large- A*' limit (long dashed). The 
spectra are COBE normalized and the cosmological parame- 
ters are S7a — , h — 0.5. 



tions of other cosmological parameters, like the addition 
of massive neutrinos, hot dark matter, which is not con- 
sidered here, will result in a change of the spectrum on 
small scales but will not resolve the discrepancy at large 
scales. 

Nonetheless, scale dependent biasing may exist and 
lead to a non-trivial relation between the calculated dark 
matter power spectrum and the observed galaxy power 
spectrum. We are thus very reluctant to rule out the 
model by comparing two in principle different things, the 
relation of which is far from understood. Therefore we 
would prefer to reject the models on the basis of pecu- 
liar velocity data, which is more difficult to measure but 
most certainly not biased. 



C. Bulk velocities 

To get a better handle on the missing power on 20 to 
100/i~^Mpc, we investigate the velocity power spectrum 
which is not plagued by biasing problems. The assump- 
tion that galaxies are fair tracers of the velocity field 
seems to us much better justified, than to assume that 
they are fair tracers of the mass density. We therefore 
test our models against peculiar velocity data. We use 
the data by Ref. |63| which gives the bulk flow 



<jUR) 



27r2 



P{k)W{kR)dk 



(122) 



in spheres of radii i? = 10 to 60ft. -'^Mpc. These data 
are derived after reconstructing the 3— dimensional ve- 
locity field with the POTENT method (see H and ref- 
erences therein). 

As we can see from Table 0, the COBE normalized 
texture model predicts too low velocities on large scales 
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FIG. 22. Matter Power spectrum: comparison between 
data and theory. References are in the text. Data set courtesy 
of M. S. Vogeley H. 
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when compared with POTENT results. Recent measure- 
ments of the bulk flow lead to somewhat lower estimates 
like ay{R) - (230 ± 90) at i? = eO/i-^Mpc ( H), but 
still a discrepancy of about a factor of 2 in the best case 
remains. 

Including a cosmological constant helps at large scales, 
but decreases the velocities on small scales. 

If the observational bulk velocity data is indeed reliable 
(there are some doubts about this |66|), all global 0{N) 
models are ruled out. 



V. CONCLUSIONS 

We have developed a self contained formalism to de- 
termine CMB anisotropics and other power spectra for 
models with causal scaling seeds. We have applied it to 
global 0{N) models which contain global monopoles and 
texture. Our main results can be summarized as follows: 

• Global 0{N) models predict a flat spectrum 
(Harrison-Zeldovich) of CMB anisotropics on large 
scales which is in good agreement with the COBE 
results. Models with vanishing cosmological con- 
stant and a large value of the Hubble parameter 
give as ~ 0.4 to 0.5 which is reasonable. 

• Independent of cosmological parameters, these 
models do not exhibit pronounced acoustic peaks 
in the CMB power spectrum. 

• The dark matter power spectrum from global 0{N) 
models with Ha = has reasonable amplitude 
but does not agree in its shape with the galaxy 
power spectrum, especially on very large scales 
> 20h-^Mpc. 

• Models with considerable cosmological constant 
agree relatively well with the shape of the galaxy 
power spectrum, but need very high bias 6 ~ 4 — 6 
even with respect to IRAS galaxies. 

• The large scale bulk velocities are by a factor of 
about 3 to 5 smaller than the value inferred from 
11. 

In view of the still considerable errors in the CMB data 
(see Fig. |2^), and the biasing problem for the dark mat- 
ter power spectrum, we consider the last argument as 
the most convincing one to rule out global 0{N) mod- 
els. Even if velocity data is still quite uncertain, ob- 
servations generally agree that bulk velocities on the 
scale of 50h~^Mpc are substantially larger than the (50 
- 70)km/s obtained in texture models. 

However, all our constraints have been obtained as- 
suming Gaussian statistics. We know that global defect 
models are non-Gaussian, but we have not investigated 
how severely this influences the above conclusions. Such 
a study, which we plan for the future, requires detailed 



maps of fluctuations, the resolution of which is always 
limited by computational resources. Generically we can 
just say that non-Gaussianity can only weaken the above 
constraints. 

Our results naturally lead to the question whether 
all scaling seed models are ruled out by present data. 
The main problem of the 0{N) model is the missing 
power at intermediate scales, i ^ 300 — 500 or i? ~ 
(20 — 100)ft.~^Mpc. We have briefly investigated whether 
this problem can be mitigated in a scaling seed model 
without vector and tensor perturbations. In this case, 
also scalar anisotropic stresses are reduced by causality 
requirements (see Ref. [Q), and the compensation mech- 
anism mentioned in Section II is effective. For simplicity, 
we analyze a model with purely scalar perturbations and 
no anisotropic stresses at all, /tt = 0. The seed func- 
tion <I>s is taken from the texture model (numerical sim- 
ulations) and we set = —^s- The resulting CMB 
anisotropy spectrum is shown in Fig. top panel. A 
smeared out acoustic peak with an amplitude of about 
2.2 does indeed appear in this model. This is mainly due 
to the fact that fluctuations on large scales are smaller 
in this model, as is also evident from the higher value of 
e = (2.2 ± 0.2) X 10~^. But also here, the dark matter 
density fluctuations and bulk velocities are substantially 
lower than observed galaxy density fluctuations or the 
POTENT bulk flows. 

Clearly, this simple example is not sufficient and a 
more thorough analysis of generic scaling seed models 
is presently under investigation. So far it is just clear 
that contributions from vector and tensor perturbations 
are severely restricted. 
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Saskatoon 


17 


1990 


950 


630 


0.0037 








0.79 


Saskatoon 


18 


4490 


1690 


1360 


0.0037 








3.83 


Saskatoon 


19 


6930 


2770 


2140 


0.0037 




4S 




4.60 


Saskatoon 


20 


6980 


3030 


2310 


0.0037 








4.01 


Saskatoon 


21 


4730 


3380 


3190 


0.0037 




4S 




1.32 


CATl 


22 


934 


403 


232 


0.0001 








1.36 


CAT2 


23 


577 


416 


238 


0.0001 




49 




0.62 



TABLE II. The CMB anisotropy detections used in our analysis. The 3., 4. and 5. column denote the value of the anisotropy 



and the upper and lower l-cr errors respectively. The references are: Tegmark and Hamilton 1997 



|9|; Cheng et al. 1994 
Gundersen et al. 1993 | 



Cheng et al. 1996 [|T|[^:]heng et al. 1997 
Dragovan et al. 19&3~^; Masi et al 1996 



Tanaka et al. 1996 
Netterfield et al. 1996 




de Bernardis et al. 1994 
Gutierrez et al. 1997 JM; 
Scott et al. 1997 lEgT 
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Catalog 


h 


17a 


Best fit bias b 




Data points 


CfA2-SSRS2 101 Mpc 


0.5 


0.0 


3.4 


29 


24 


CfA2-SSRS2 101 Mpc 


0.8 


0.0 


2.0 


40 


24 


CfA2-SSRS2 101 Mpc 


1.0 


0.0 


1.9 


44 


24 


CfA2-SSRS2 101 Mpc 


0.5 


0.4 


3.9 


17 


24 


CfA2-SSRS2 101 Mpc 


0.5 


0.8 


9.5 


4 


24 


CfA2-SSRS2 130 Mpc 


0.5 


0.0 


5.3 


8 


19 


CfA2-SSRS2 130 Mpc 


0.8 


0.0 


3.4 


15 


19 


CfA2-SSRS2 130 Mpc 


1.0 


0.0 


3.4 


16 


19 


CfA2-SSRS2 130 Mpc 


0.5 


0.4 


5.6 


5 


19 


CfA2-SSRS2 130 Mpc 


0.5 


0.8 


11.1 


4 


19 


LCRS 


0.5 


0.0 


3.0 


71 


19 


LCRS 


0.8 


0.0 


1.8 


96 


19 


LCRS 


1.0 


0.0 


1.6 


108 


19 


LCRS 


0.5 


0.4 


3.7 


33 


19 


LCRS 


0.5 


0.8 


8.7 


40 


19 


IRAS 


0.5 


0.0 


2.3 


102 


11 


IRAS 


0.8 


0.0 


1.3 


131 


11 


IRAS 


1.0 


0.0 


1.3 


140 


11 


IRAS 


0.5 


0.4 


2.8 


70 


11 


IRAS 


0.5 


0.8 


6.3 


9 


11 


IRAS 1.2 Jy 


0.5 


0.0 


4.2 


56 


29 


IRAS 1.2 Jy 


0.8 


0.0 


2.9 


92 


29 


IRAS 1.2 Jy 


1.0 


0.0 


2.9 


99 


29 


IRAS 1.2 Jy 


0.5 


0.4 


4.3 


39 


29 


IRAS 1.2 Jy 


0.5 


0.8 


6.7 


28 


29 


APM 


0.5 


0.0 


3.3 


1350 


29 


APM 


0.8 


0.0 


1.8 


1500 


29 


APM 


1.0 


0.0 


1.7 


1466 


29 


APM 


0.5 


0.4 


3.5 


1461 


29 


APM 


0.5 


0.8 


6.2 


1500 


29 


QDOT 


0.5 


0.0 


4.3 


32 


19 


QDOT 


0.8 


0.0 


2.9 


44 


19 


QDOT 


1.0 


0.0 


2.9 


46 


19 


QDOT 


0.5 


0.4 


4.3 


25 


19 


QDOT 


0.5 


0.8 


7.3 


14 


19 



TABLE III. Analysis of the matter power spectrum. In the 
model parameters. In cols. 4 and 5 we give the bias parameter 
shows the number of 'independent' data points assumed in the 



first column the catalog is indicated. Cols. 2 and 3 specify the 
inferred by minimization as well as the value of x^- Col. 6 
analysis. 



R 


<yv (R) 


A„ 


h ^ 0.5 


h = 1.0 


Ha = 0.8 


10 


494 


170 


145 


205 


86 


20 


475 


160 


100 


134 


78 


30 


413 


150 


80 


98 


70 


40 


369 


150 


67 


78 


65 


50 


325 


140 


57 


65 


61 


60 


300 


140 


50 


56 


57 



TABLE IV. Bulk velocities: Observational data from [g3| and theoretical predictions. Ay estimates the observational 
uncertainty. The uncertainties on the theoretical predictions are around ~ 30%. The models SIa ~ with h — 0.5 and h — 1 
as well as = 0.8, h = 0.5 are investigated. 
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APPENDIX A: COMPLETE DEFINITIONS OF 
GAUGE-INVARIANT PERTURBATION 
VARIABLES 

In this Appendix we give precise definitions of all the 
gauge-invariant perturbation variables used in this pa- 
per. These definitions, their geometrical interpretation 
and a short derivation of the perturbation equations can 
be found in |^,|l3|. We restrict the analysis to the spa- 
tially flat case, K — 0. We define the perturbed metric 

by 



g = g + a'^h , 



(Al) 



where g denotes the standard Friedmann background, a 
is the scale factor and h denotes the metric perturbation. 



1. Scalar perturbations 

Scalar perturbations of the metric are of the form 

h^s) ^ ^2A{dtf + 2iBkjdtdx^ + 2{Hl + -HT)5^jdx'^dx^ 

3 

-2k-^HThkjdx'dx^ . (A2) 

Computing the perturbation of the Ricci curvature scalar 
and the shear of the equal time slices, we obtain 



1 



6R = 4a~^fc^7^ , with TZ = Hl + -Ht , (A3) 



with 



cr = k'^Hr - B 



The Bardeen potentials are the combinations 



$ = 7^ - (d/a)fc-V 

* = A - k^'^[{a/a)(j - a] 



(A4) 
(A5) 



(A6) 
(A7) 



They are invariant under infinitesimal coordinate trans- 
formations (gauge transformations). 

To define perturbations of the most general energy mo- 
mentum tensor, we introduce the energy density p and 
the energy flux u as the time-like eigenvalue and normal- 
ized eigenvector of Tj^, 

T/u^' = -pu" , = -1 . 

We then define the perturbations in the energy density 
and energy velocity field by 



P = -p{l + 5) 



(A8) 



(A9) 



vP is fixed by the normalization condition, = a~^{\ — 
A). In the 3-space orthogonal to u we define the stress 
tensor by 



t"" EE P^' ^P" pT'^'^ 



(AlO) 



where P ~ u®u + g is the projection onto the sub-space 
of TAi normal to u. It is 



r° - r° - r' - n 



The perturbations of pressure and anisotropic stresses 
can be parameterized by 

r/' = p[{l + ttl)5P + TT, , with ttI^O. (All) 

For scalar perturbations we set 



u 
and 



nf ' = i-k-^k'kj + -s'^)n . 



Studying the behavior of these variables under gauge 
transformations, one finds that the anisotropic stress po- 
tential n is gauge invariant. A gauge invariant velocity 
variable is the shear of the velocity field. 



(Sm) 



1 



{k^^hkj Sijtt'^V , with V = v- k^^Hj 



(A12) 

There are several different useful choices of gauge in- 
variant density perturbation variables, 

Ds=5 + 2,{l + w){a/a)k-^(j (A13) 
Dg=5 + Z{l + w)n = D, + i{l + w)<^ (A14) 
D = Ds + -i{l + w){d/a)k-^V . (A15) 

In this work we mainly use Dg. Here w — p/p denotes 
the enthalpy. Clearly, these matter variables can be de- 
fined for each matter component separately. For ideal 
fiuids like CDM or the baryon photon fluid long before de- 
coupling, anisotropic stresses vanish and ttl = {c1/w)5, 
where Cs is the adiabatic sound speed. 

Also scalar perturbations of the photon brightness, l^^^ 
are not gauge invariant. It has been show that the 
combination 



,is) 



Aik ^nPkjCj 



(A16) 



is gauge invariant. This is the variable which we use 
here. In other work | |67[ | the gauge invariant variable 
<d = M + ^ has been used. Since $ is independent of the 
photon direction n this difference in the definition shows 
up only in the monopole, Cq. But clearly, as can be seen 
from Eq. (A16), also the dipole, Ci, is gauge dependent. 

The brightness perturbation of the neutrinos is defined 
the same way and will not be repeated here. 
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2. Vector perturbations 

Vector perturbations of the metric are of the form 

h^^'> = 2Bjdx'dt + ik^^ikiHj + kjHi)dx^dx^ , (A17) 

where B and H are transverse vector fields. The simplest 
gauge invariant variable describing the two vectorial de- 
grees of freedom of metric perturbations is S, 



— fc ^Hj B-j 



(A18) 



Vectorial anisotropic stresses are gauge invariant. They 
are of the form 

tt\J^ = ik-\k,I{i + kiUj) . (A19) 

The vector degrees of freedom of the velocity field are 
cast in the vorticity 



ui:j - Uj-i = ia{kjUJi - kiLOj) 
with LOj = Vj — Bj . 



(A20) 
(A21) 



Vector perturbations of the photon brightness are gauge- 
invariant. To maintain a consistent notation, we denote 
them by M^^'> . 



3. Tensor perturbations 



We define tensor perturbations of the metric by 



/i(^) = 2H,jdx'dx^ 



(A22) 



where Hij is a traceless transverse tensor field. 

The only tensor perturbations of the energy momen- 
tum tensor are anisotropic stresses, 



Tensor perturbations of the photon brightness are de- 
noted M^'^l 

Clearly, all tensor perturbations are gauge-invariant 
(there are no tensor type gauge transformations). 



APPENDIX B: THE CMB ANISOTROPY 
POWER SPECTRUM 



Here we derive in some detail Eqs. (|75|), ( p4\ ) and (91). 

CMB anisotropics are conveniently expanded in spher- 
ical harmonics: 6T{n)/To = J2im^im.Ym(.^)- The coeffi- 
cients aim are random variables with zero mean and ro- 
tationally invariant variances, = (| aim P)- The mean 
(over the ensemble) correlation function of the anisotropy 
pattern has the standard expression: 

|^(ni)^(n2)) ^ ^ Y.(^e + l)CePiicose) (Bl) 



where cos = ni • n2. To find Eq. 
transform normalization 



we use the Fourier 



/(k) 



1 
V 



/(x) exp(ik • x)d3 



(B2) 



with some normalization volume V. Assuming that en- 
semble average can be replaced by volume average then 
implies 

/6T. .ST. ,\ 1 f 3 ST. .ST. . 

= 1^ I d'fc^(k,ni)^(k,n2) . (B3) 
[2tiY J To To 

Inserting our ansatz (^J) for ^ = jM., and using the 
addition theorem for spherical harmonics, we have 



ST ST 
77r(ni)77r(n2) 
Jo -to 



1 ^ (-l)(^-^')r,™(ni)yA„,(n2)x 



edkdn^Y.l.^ (k)r,,„, (k) {cJicrDik) 



-^5](2^+ l)P,(ni • n2) / edk{a,al){k) , (B4) 



from which we can read of Eq. (75) 



For vector and tensor fluctuations, the ansatz (77) and 



( p6| ) must be taken into account. With the same ma- 
nipulations as above the correlation function of CMB 
anisotropics induced by vector modes reads: 



ST ST \ 
-to Jo / 

(A23) /'d3fc^n,i,£2(k,n'i,n'2)Pn(A*'i)T£2(A^'2) (B5) 



where the primes indicate that the quantity is considered 
in the reference system where k is parallel to the z axis 
and ( !§], @) 



[(^l/l'^l/2 >COs(0i)( 



+ (4!^i'^2!f2*>sin('/''i)sin((/)2)] • 
Assuming statistical isotropy which implies 

iWZ^?) = {WZ^?) and 
we obtain 



(B6) 



= , 
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' 6T jn,) ST in2) 
V To To 

^^(2^i + l)(2£2 + l)H)("- 



128773 



(B7) 



eu2 



J Tin\,n',){a^^,y^>)Pn{t^\)Pait^',)d'k , (B8) 



where 



T = ni • n2 - AiiM2 



(B9) 



Using the recursion formula for Legendre polynomials 
and the addition theorem for spherical harmonics, we 
find after some manipulations 



C. 



Btt 



k^dk- 



(2^+1)2 



(BIO) 



For the correlation function of the CMB anisotropies 
from tensor modes our ansatz (Bq) gives 



ST ST 
77r(ni)77r(n2) ) = 
-to Jo 



1287r 



with 



na/2 - i-if-'\2£l + 1)(2£2 + 1)(1 - (^^'1)^) x 

(1 - if^'^fWZlKjZn) cos(20;) cos(20'2) + 



+ (a(.^]ia;^]2>sin(2^'i)sin(20'2)]. 
Here, statistical isotropy leads to 
'dT{ni)ST{n2)' 



(B12) 



To To 
^^(24 + l)(24 + l)H)("- 



1287r3 



ai2 



J TK,n'2)(aJ']iaJ'i2*)FnK)P.2(A*'2)rf'fc (B13) 
where 

T = [2(ni . n2 - /.Vi)^ - (1 - {f^m • 

(B14) 

With straightforward but somewhat cumbersome manip- 
ulations, applying the recursion formula for Legendre 
polynomials and the addition theorem for spherical har- 
monics, we then obtain 

- 8n{e-2)\i (2£+l)2 



with 



(T) 



2(2i?+l)a^7 



JT) 



2£~ 1 (2£- l)(2£-f 3) 2e + 3 



(B16) 



The formulas (pl),(^10|) and ( |B15|) are used to determine 
the CMB anisotropy spectrum. 
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